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Abstract 

In this letter we show that the late-time scaling state in spinodal decomposition is not unique. 
We performed lattice Boltzmann simulations of the phase-ordering of a 50%-50% binary mixture 
using as initial conditions for the phase-ordering both a symmetric morphology that was created by 
symmetric spinodal decomposition and a morphology of one phase dispersed in the other, created 
by viscoelastic spinodal decomposition. We found two different growth laws at late times, although 
both simulations differ only in the early time dynamics. The new scaling state consists of dispersed 
droplets. The growth law associated with this scaling state is consistent with a L ~ t 1 ^ 2 scaling law. 



I. INTRODUCTION 



II. NUMERICAL METHOD 



The phase-ordering in binary fluids following spinodal 
decomposition has beenJjhe focus of considerable inter- 
est in the recent pastErfcl. A number of different scal- 
ing regimes with different growth laws have been iden- 
tified with the help of scaling arguments, but little em- 
phasis has been placed on morphological effectsoS. Un- 
usual asymmetric morphologies are often seen WrUie early 
stages of viscoelastic spinodal decompositiomTla, but it 
has been argued that viscoelastic effects are unimpor- 
tant at late times and therefore only one true late-time 
phase-ordering scaling state existsB. The term "late-time 
scaling state" is not clearly defined since a number of scal- 
ing states that cross over at certain length-scales have 
been predicted and observed. In the present context, we 
take a "late-stage" scaling state to mean a fully phase- 
separated state in which the time dependence of a typical 
length scale L(t) has converged to a growth law of the 
form L(t) ~ t a . Even though the existence of these, scal- 
ing states in a strict sense has been questioned™, the 
concept remains useful to understanding phase-ordering 
dynamics in binary fluids. 

The argument of the unimportance of viscoelasticity 
in the late- stage growth, as well as many of the scaling 
arguments, rely on an assumption that is often implic- 
itly made: for the phase-ordering of symmetric binary 
fluid mixtures only one universal late-time scaling state 
exist sa. In this letter, however, we show that this as- 
sumption is not valid, and that therefore the early time 
dynamics are very important in selecting the late-time 
scaling state. The possibility of the non-uniqueness of 
the scaling state was first suggested by A. Rutenberg, 
although at tfi£_time no persuasive numerical evidence 
could be foundE£l. 



For the simulations we use the viscoelastic two- 
component lattice Boltzmann simulation introduced in 
an earlier papertil. Briefly, in lattice Boltzmann simula- 
tions densities fa that are associated with velocities Vi are 
streamed on a lattice according to the lattice Boltzmann 
equation 



/i(x + v*AM + At) = 



(1) 



where ff is an equilibrium distribution, A^ is a colli- 
sion matrix, and v^At is a lattice vector. The veloc- 
ity set for our simulation consists of 17 velocities given 
by {(0,0), (1,0), (0,1), (-1,0), (0,-1), (1,1), (-1,1), 
(-1,-1), (i,-i), (1,0), (0,1), (-1,0), (0,-1), (1,1), 
(-1,1), (-1,-1), (1,-1)}. Note that the last 8 veloci- 
ties are the same as the previous eight velocities. This 
duplicity allows the simulation to have two independent 
stresses which represent a viscoelastic and a purely vis- 
cous contribution to the total stress tensor. The two 
contributions are used to produce a Jeffrey's model for 
the stress (see eqn. (||)). The algorithm is required to 
conserve mass and momentum, but not energy. Energy 
conservation is replaced by a condition of constant tem- 
perature. The macroscopic density p and velocity u are 
defined as 



(2) 



To simulate a two-component mixture we introduce a 
second lattice Boltzmann equation for densities gi 



gi (x + ViAt,t + At) = 

g(x,t) + — [^(x,t)-</i(x,t) 



(3) 



1 



where we choose a diagonal collision matrix with a single 
relaxation time r. These densities are only defined on the 
first nine velocities v^. The density difference (j) = p\— p2 
of the two components is given by 



9i 



(4) 



and the total density is given by p = p\ + p2 . By choosing 
appropriate equilibrium distributions and an appropriate 
collision matrix, we ensure that the following partial dif- 
ferential equations are being simulated up to second order 
in the derivatives but assuming that the relaxation of the 
viscoelastic stress a is slow (0 ~ l/v^) : 

d t p + d x (pu) =0 (5) 

pd t u + pu.Vu = -<9 X P + d x (a v + a) (6) 
cr v = ^oo(V(pu) + (V(pu)) T - V.utf) 

+^ooV.U(5 (7) 

a + 0(0)<t (1) = -(i/ o (0) - ^oo)(V(pu) + (V(pu)) T ) (8) 

d t <j) + 9 x (0u) = DV 2 /i + V.((0/p)V.(P - a)) (9) 

where (5 is the identity matrix, a is the viscoelastic stress 
tensor, am = c^a + u.Va — cr.(Vu) — (Vu) t <j is its 
upper convected derivative, P = 0.5p + 0.007(V(/>V(/> — 
O.5V0.V05 — (j)V 2 (j)5 is the pressure tensor, and p = 
-0.550/p + 0.25 ln((p + 0)/(p - 0)) - 0, OO7V 2 is the 
chemical potential. The parameters v^, £oo, and are 
determined by the eigenvalues of the collision matrix. 
The values for the parameters were At = 1, r = 1, 
D 0.5, 0.31, and 0.01. For the low- 

viscosity phase we used = 0.055, v§ = 0.013 and for 
the viscoelastic phase = 39.5 and v$ = 1.97. For the 
symmetric simulations of Figure [| we used = 0.055 and 
V( } « z/oc = 0.075. 



III. SIMULATIONS 

We performed simulations of critical spinodal decom- 
position of a viscoelastic binary mixture in two dimen- 
sions where one component is much more viscoelastic 
than the other. These simulations lead to the usual mor- 
phologies in which the viscoelastic phase-is connected and 
the less viscoelastic phase is dispersed!^. We performed 
our simulations on a 256 2 lattice and after about 1000 
iterations the less viscoelastic phase is completely dis- 
persed, although the domains are still highly deformed. 
We used this morphology as an initial condition for a sim- 
ulation where we make both components purely viscous 
to examine the effects of initial conditions that are not 
symmetric on symmetric binary fluid mixtures. This al- 
lows us to distinguish the effect that the morphology cre- 
ated by viscoelastic phase separation has from the effect 
of viscoelasticity itself in the late-time phase-ordering 
process. 
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FIG. 1. A comparison of the phase-ordering of two identi- 
cal symmetric binary mixtures after different early time spin- 
odal decomposition. In (a) a phase-ordering is seen from 
an initial morphology generated by viscoelastic phase separa- 
tion in which the light component is dispersed (the originally 
low- viscous component). In (b) we see the usual symmetric 
phase separation. Fig.(c) shows the length scales L(t) for (a) 
L~t 1/2 and (b) L~t 2/3 . 
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It is usual to assume that the phase-ordering process 
is independent of the details of the early time spinodal 
decomposition, although this assertion has, to our knowl- 
edge, only been tested with respect to different random 
initial conditions. If this assumption were true we would 
expect the domain growth to lead to a reconnection of 
the dispersed domains that would lead to a morphology 
and a growth law that are identical to a system which 
started with the same viscosities for both components. 

In Figure [j] we see a comparison of the phase-ordering 
from an initially dispersed phase (a) and a symmetric 
initial condition (b). The morphologies for both sys- 
tems are shown after 1000, 2000, 4000, and 8000 iter- 
ations. We see that the phase-ordering of the morphol- 
ogy with a dispersed phase leads to an even more pro- 
nounced dispersed phase where almost all domains are 
circular at late times. We see that droplet coalescence 
occurs frequently with only few domains vanishing due 
to the evaporation-condensation mechanism underlying 
Oswald ripening. The droplet coalescence, however, is 
not frequent enough to change the connectivity of the do- 
mains. Instead, we observe that domains become more 
circular on average, suggesting that even for very long 
times we do not expect a transition to a bi-continuous 
morphology. In Figure |l|(c) we see that the growth law 
for the droplet phase is L(t) ~ t 1 ! 2 and is smaller than 
the L(t) ~ t 2 / 3 seen for the symmetric phase-ordering 
shown in Figure 0(b), but also faster than the L(t) ~ t 1 ! 3 
usually seen in the coarsening of off-critical droplet mor- 
phologies of symmetric binary mixtures. (The data for 
the well established scaling law of L ~ t 2 ! 3 for the sym- 
metric phase-ordering in Figure p](c) is more noisy since 
it relies on a single simulation on a 256 2 lattice where as 
the data for the droplet morphology has been obtained 
by simulation on a 1024 2 lattice.) We suggest that the 
high density of droplets enhances the droplet coalescence 
an leads to a faster growth of L ~ t 1 ! 2 . We measure L(t) 
as the inverse of the length of the interface. We should 
also point out that a length scale derived as L ~ 1/ y/N, 
where N is the number of domains, scales as L ~ til 2 as 
well, suggesting that it may be a true scaling state! 2 ]. 

The existence of this second scaling state, distinct from 
the bi-continuous scaling state, is important to under- 
standing the late-time regime of viscoelastic phase sep- 
aration because, even in the absence of viscoelastic ef- 
fects, we observe a droplet-morphology evolving from the 
initial morphology created by viscoelastic phase separa- 
tion. This result is also important for practical reasons in 
processes where late-time morphologies need to be con- 
trolled. It is well known that mixing of high- viscosity and 
low- viscosity components by means of mechanical agita- 
tion leads to morphologies where the high- viscosity phase 
is dispersed in dropletsll3 for volume fractions of the low- 
viscosity component of much less than 50%. This-effect is 
enhanced if the high- viscous phase is viscoelastidi-l This 
leads us to consider what the late-time morphology of a 
phase-ordering system with an early morphology created 
by mechanical mixing would be. 




-9.5 [/ i i i i i i i i i : 

7 8 9 10 11 

ln(t) 

(c) 

FIG. 2. The phase-ordering of a binary mixture of a vis- 
coelastic phase (dark color) and a low viscosity Newtonian 
phase (light color) is shown in (a) and (b). (a) corresponds to 
a usual viscoelastic phase separation of a 50%-50% mixture 
and (b) shows the phase-ordering after the viscoelastic and 
low- viscosity domains have been interchanged. Morphologies 
are shown after 2000, 4000, 8000, and 16000 iterations. The 
scaling behavior of the two mixtures is shown in (c) where for 
(a) L ~ t 2 and for (b) L ~ 0.4. 
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In order to answer this question we performed a simu- 
lation of viscoelastic phase separation and after a droplet 
morphology had been formed at 1000 time-steps, we in- 
verted the properties of the two components. We used 
this state as a model for the morphologies of dispersed 
droplets of the viscoelastic phase in a matrix of the low 
viscosity Newtonian phase that is typical of mechani- 
cally mixed morphologies. We then continued the simu- 
lations and observed the phase-ordering behavior of the 
new morphology. The results of these simulations are 
shown in Figure ||. In Figure |^(a) the morphology for 
a viscoelastic phase separation of a 50%-50% mixture is 
shown. We should emphasize that the phase-ordering 
(see eqn. (^|)) does not have a composition-dependent 
diffusion constant and therefore no domain shrinkage is 
observed in these simulations. Domain shrinkage can be 
a very slow process in viscoelastic phase separation that 
prolongs the spinodal decomposition process and makes 
it more difficult to differentiate the early-stage decompo- 
sition and the late-time phase-ordering processes. This 
does not reduce the validity of our results, however, since 
we are only interested in the late-time behavior when the 
domain shrinkage is completed. 

The morphologies shown in Figure ^|(b) are of a simu- 
lation where at 1000 iterations, after the dispersion was 
achieved, the viscoelastic and the low viscosity compo- 
nents were exchanged. We see that the morphology re- 
mains a droplet morphology, albeit now of the viscoelas- 
tic phase. Comparing Figures | (a) and (b) we see 
that the most important factor in determining the fi- 
nal morphology is the early stage dynamics, and that 
phase-ordering scaling morphologies of both dispersed 
viscoelastic and dispersed low viscosity domains exist. 

From Figure ||(c) we see that each scaling state has 
a different growth law. The morphology of dispersed 
viscoelastic domains grows as L ~ t 0A whereas the 
morphology of dispersed low viscosity domains grows as 
~ t 0,2 after near circular droplets have been formed. An 
anomalously slow growth in viscoelastic phase separation 
has first been observed experimentally by Tanakat^l who 
found L ~ t 015 for a system of high molecular- weight 
polystyren/ diethyl malonate (4.0 wt.%). 

These simulations also emphasize the difference be- 
tween a morphology after spinodal decomposition and 
after mechanical mixing. The morphology after spin- 
odal decomposition is a dispersed low viscosity phase, 
whereas the state after mechanical mixing has a dispersed 
viscoelastic phase. Subsequent phase-ordering does not 
change the connectivity of these states, in agreement with 
the conventional wisdom that there is a profound differ- 
ence in states produced by spinodal decomposition and 
mechanical mixing. 



IV. CONCLUSIONS 

In this letter we have shown that more than one scal- 
ing state exists for late-time spinodal decomposition and 
that the early time behavior of a phase-separating binary 
mixture is very important in selecting one of these scal- 
ing states. We have also explained the difference between 
viscoelastic phase-ordering states after spinodal decom- 
position and mechanical mixing. Our results show that 
the volume fraction and the physical properties of a mix- 
ture do not select a morphology by themselves, but that 
the morphology of the initial state is of paramount im- 
portance. This is why viscoelastic phase separation can 
lead to unusual late-time scaling states even when vis- 
coelasticity is no longer important at large length scales. 
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